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Motivated by recent "circuit QED" experiments we investigate the noise properties of coherently 
driven nonlinear resonators. By using Josephson junctions in superconducting circuits, strong non- 
linearities can be engineered, which lead to the appearance of pronounced effects already for a low 
number of photons in the resonator. Based on a master equation approach we determine the emission 
spectrum and observe for typical circuit QED parameters, in addition to the primary Raman-type 
peaks, second-order peaks. These peaks describe higher harmonics in the slow noise-induced fluc- 
tuations of the oscillation amplitude of the resonator and provide a clear signature of the nonlinear 
nature of the system. 
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I. INTRODUCTION 

Several recent circuit QED experiments based on 
Josephson junctions in superconducting circuits were 
concerned with the properties of driven nonlinear res- 
onators |l|. Josephson junctions allow engineering strong 
nonlincaritics, which can lead to pronounced quantum ef- 
fects, such as non-classical photon number distributions 
k|. At the same time, without leaving the quantum 
regime, the damping of the resonator can be made suf- 
ficiently strong to measure the radiation emitted by the 
resonator. State-of-the-art measurements give access to 
the spectral density of the emitted radiation S(ui). This 
yields detailed information about the quantum and clas- 
sical fluctuations of the resonator [J-Q ■ 

Some of the experiments with coherently driven res- 
onators based on a Josephson junction were concerned 
with the development of a Josephson bifurcation am- 
plifier to be used as a high-contrast readout device for 
superconducting qubits j6j. The Josephson bifurcation 
amplifier makes use of the dynamical bistability induced 
by a linear periodic driving with frequency ujp close to 
the oscillator eigenfrequency ujq. In the readout process, 
the two qubit states are mapped onto the two stable vi- 
brational states of the resonator, which can be easily dis- 
tinguished, since they differ strongly in amplitude and 
phase. The measurment backaction depends on the os- 
cillator power spectrum p], @] . 

The stable vibrational states are described by nearly 
sinusoidal oscillations with frequency u)p. They are sep- 
arated by a dynamical barrier in phase space [9|. The 
amplitude and phase of the forced oscillations display 
quantum and thermal fluctuations. For typical parame- 
ters, the timescales for these fluctuations are set by the 
detuning 5u> = ujp — ujq between the coherent driving and 
the resonator and by the relaxation rate T. 

Here we consider the regime, where the fluctuations 
are small on average and the oscillator stays in the vicin- 



ity of one of the stable solutions for times much longer 
than the fluctuations characteristic timescales, T _1 and 
\6u>\ . Many theoretical investigations have focused on 
rare large fluctuations which allow overcoming the dy- 
namical barrier and induce the switching between the 
stable solutions [9, [l0| or on tunneling [IlJ. The deep 
quantum regime, where the average number of excita- 
tions in the resonator is small and quantum fluctuations 
are large has also attracted considerable attention jT3 |. 
The power spectrum in this regime has been investigated 

in m. 

Although the emission spectrum S(w) in principle car- 
ries information about rare large fluctuations [5|, Il4j . 
most of the radiation is emitted as a consequence of 
small fluctuations. Therefore, the standard approach for 
computing the spectrum S(uA is to linearize the motion 
around the stable solutions |{J 0, ■ The small fluc- 
tuations display damped harmonic oscillations with fre- 
quency ~ \6u\. The decay of the fluctuation correlations 
is accompanied by the emission of radiation with frequen- 
cies u>f — Slo and ujp + 5uj, since the slow oscillations 
are superimposed to the fast oscillations with the driving 
frequency. The resulting emission spectrum is reminis- 
cent of the Raman spectrum of a diatomic molecule: be- 
side the main Rayleigh peak at frequency ujp, which de- 
scribes radiation emitted coherently, it has two Raman- 
type side-peaks with Raman shift ±|<5u;|. In a quantum 
mechanical picture, the quasi-energy level spacing, i.e., 
the distance between the eigenstates of the Hamiltonian 
in a frame rotating with the driving frequency, is ap- 
proximatively given by /i|<5u;|, and the Raman lines are 
induced by transition between nearest-neighbor quasi- 
energy eigenstates. In this context, the same driving field 
defines and probes the quasi-energy spectrum. 

In this paper, we are interested in a regime where typ- 
ical fluctuations are intermediate in size: they are too 
weak to induce switching between the two solutions, so 
that the resonator stays locked to one of the stable so- 



lutions for a long time, but they are strong enough to 
lead to nonlinear dynamics, so that one has to go beyond 
the simple Raman picture for the emission spectrum. It 
has been shown that quantum fluctuations of intermedi- 
ate size may give rise to a fine structure in the spectrum 

In this work we point to two additional features of the 
emission spectrum, which are consequences of the nonlin- 
ear nature of the quasi-energy Hamiltonian and should be 
observable in experiments [16| : i) Intermediate strength 
nonlinear fluctuations are not simply damped sinusoidal 
oscillations but have higher harmonics yielding additional 
peaks in the resonator emission spectrum at frequencies 
lof ± n5uj, n = 2, 3, • • • . These higher order peaks resem- 
bles the peaks in the power spectrum of a static oscillator 
with a small cubic anharmonicity close to integer multi- 
ples of its eigenfrequency [17| . At the quantum level, the 
higher harmonics derive from transitions between quasi- 
energy levels which are not nearest neighbors, ii) More- 
over, fluctuations in an asymmetric quasi-energy poten- 
tial yield an additional broad peak at the frequency lof- 
Unlike the fine structure of the emission spectrum, the 
effects wc consider here are not restricted to the quantum 
regime but extend to the classical regime where thermal 
fluctuations are dominant. 

The paper is organized as follow: In Sec. [TT] wc intro- 
duce the model for the driven nonlinear resonator and 
the approach based on a master equation. In Sec. lIIII we 
review the mean field and small fluctuations theory for 
the oscillator. In Sec. II VI we outline our analytical calcu- 
lation of the emission spectrum. In Sec. |V] we compare 
our analytical finding with numerical results which are 
exact within our master equation approach. 



II. THE MODEL 

We consider a Duffing oscillator with eigenfrequency 
wq, coordinate q, momentum p, and Hamiltonian 



H(t) 
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WqV 



7 4 



Fqcos(ujFt) . (1) 



We assume that the detuning Slo = cop — <-0q of the 
driving is small, \5lo\ <C ujq, and that the nonlinearity 
for typical values of q satisfies the condition |7|(g 2 ) <C 
Wq. In this regime, the oscillator displays fast sinusoidal 
oscillations, (q(t)) sa A(t) cos[w pt + <p(i)], with amplitude 
A(i) and phase 4>{t) which vary slowly on the timescale 
\Su\-\ 

It is convenient to study the oscillator dynamics in a 
frame rotating with the driving frequency. I.e., we per- 
form a unitary transformation H = U^HU — iHU'U with 
U(t) = exp[— iuipa^at], where a = (u>oq + ip) /y/2fru>o- In 
the rotating frame, the ladder operators a and a* vary 
slowly, (a(t)) « ^cj a 2hA(t) exp(— i<fi{t)). Wc therefore 
use the rotating wave approximation (RWA) (neglecting 
fast oscillating terms with frequencies 2ujp and 4uip) and 



arrive at the time-independent Hamiltonian 

H ps -hSu a f a H a^a^a + 1) (a + a f ) , (2) 



with K = 3h-y/2uj% and / = F/^/2HlJo~. 

We consider Markovian dissipation due to a linear cou- 
pling to a bosonic bath. When the oscillator relaxation 
rate is much smaller than its eigenfrequency, r < wo, 
and the bath spectral density close to wo is smooth, the 
dissipative dynamics of the system can be described by 
a simple master equation in the Lindblad form 

p = Cp = --[#, p] + r(l + n)V[a]p + TnV[a)]p, (3) 

where the Lindblad supcropcrator is defined through 
V[0}p = 20pO J < - O^Op - pOO\ and n = ( e h ^/k B T _ 
1) _1 is the oscillator distribution in the absence of driv- 
ing. 

For long times the resonator relaxes to a stationary 
state, satisfying Cp st = 0, from which we obtain mean 
values (O)st = Tr(p st O) in the rotating frame. In par- 
ticular, it allows us to calculate the oscillation amplitude 
A st = VW^(a)st. 

The Lindblad master equation §S§ gives also access to 
the full time-evolution of the system, and wc can use 
it to calculate correlation functions (0(t + Ai)0'(t)) = 
Tr{Oe CAt O' p(t)} [lq |, where C is the superoperator de- 
fined in Eq. Q. For long times t, the system is in the 
stationary state p st , so that the correlation function de- 
pends on the time difference At only, (0(t + Ai)O'(i)) = 
Tr{Oe £At OV st } = (0(At)0') st . 

Here, wc arc specifically interested in the emission 
spectrum S(u), i.e., the spectral density of the photons 
emitted by the driven resonator, 



S(u) = 2 Re / dt (a^ (t)a) st e- i(u - UF)t 
Jo 



(4) 



The above definition takes into account that the correla- 
tion function is computed in a frame rotating with fre- 
quency UJp. 

The emission spectrum S(lo) consists of two distinct 
contributions coming from the finite mean value (a) s t and 
from the fluctuations ((<v(t) — (a) s t)(a — (a) st )) s t of the 
operator a (which describes the oscillation amplitude). 
The latter yield broad peaks with width of the order of 
the oscillator relaxation rate T and height proportional 
to the noise intensity. On the other hand, the finite mean 
value (a) s t yields a sharp peak at the driving frequency 
Lop. For a strictly monochromatic driving, this peak has 
the shape of a delta function, Tr(a}g t 6(<j0 — u>f). In the 
following, we subtract this contribution from the emission 
spectrum S(u>) and focus only on the broad peaks due to 
the fluctuations of a. 

In Fig.[T] we show the emission spectrum S(ui) for a pa- 
rameter choice in the regime of operation of a Josephson 
bifurcation amplifier [l|. In this example, the lineshapc 
of S(u>) consists of four distinct Lorentzian peaks: two 
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Figure 1. Emission spectrum of the driven nonlinear resonator 
in the regime of high-amplitude oscillations. The parameters 
are: 5co = -75MHz, K = -3MHz, / = 270MHz, F = 3MHz 
and n — 1. The dashed lines show the second-order peaks 
scaled up by a constant factor. 



high peaks of width T for u> — ujf close to ±8lo and two 
lower peaks of width 2T for u — tup close to ~2Su> and 
0. Such peaks are consistent with weakly damped oscil- 
lations of the fluctuation of a with period Sui, which are 
mostly harmonic but have substantial second harmonic 
and average over a period 2tt/Slu. A linearized theory of 
the fluctuations around the mean-field solutions [f|, [l5( 
yields damped harmonic oscillations for the fluctuations 
of a, but does not capture the second order peaks for 
lu — top w — 2Suj, 0. In order to go beyond this limitation 
we will compute in Section IIVI the emission spectrum by 
taking into account the nonlinear nature of the fluctua- 
tions. 



III. THE MEAN FIELD SOLUTIONS AND 
SMALL FLUCTUATIONS 

In this section, in order to prepare for the calculation 
of the emission spectrum, we review the mean field and 
small fluctuation theories for the Duffing oscillator [9|, [lj, 

M. 



A. Dynamical bistability 

For concreteness, we focus on soft nonlinearities, K < 
(which apply to the Joscphson nonlincarity), and red 
detuned driving, 8uj < 0. We switch to dimensionlcss 
variables by introducing the time r = \8u>\t, the effective 
Planck constant 



K 

2Ki 



3hj 




Figure 2. 3D-plot of the function g(Q, P) for a fixed value of 
the scaled driving /3 = 0.034. The minimum corresponds to 
the small amplitude state and the maximum corresponds to 
the high amplitude state. 



and the slowly varying quadratures Q and P, 



Q = \l^(a- 



at), 



P = 



-i{a) -a) 



with canonical commutator [Q, P] 
less Schrodingcr equation reads 



(0) 
iX. The dimension- 



iXdMQ, t) = g(Q, P = - 4 A9q)V(Q, t) . (7) 
where the quasi-energy Hamiltonian 

g = -(Q 2 + P 2 -lf/A+^Q. (8) 

depends on the parameter j3 = f 2 K / (A\8uj\ z ) only and is 
time- independent. However, it can not be written as a 
sum of kinetic and potential energy. For hard nonlinear- 
ity, K > 0, and blue detuning, Slo > 0, the quasienergy 
surface g has opposite sign |j|. 

The values of quasi-energy Hamiltonian g(Q, P) as 
function of the variables Q, P are shown in the 3D plot 
in Fig. [2] for a fixed value of j3 = 0.034. It has three sta- 
tionary points corresponding to a minimum, a maximum, 
and a saddle point. These stationary points coexist for 
< j3 < A/27 . The maximum and the minimum are 
located at 



(Q, P) = (Q h = cos 0/y/3 + sin 0, P h = 0) (9) 



and 



(Q, P) = (Qi = cos 0/Vs - sin 9, P\ = 0), 



(10) 



(5) 



respectively. 

Here, the angle 8 is given by 8 = (ir — 
arctan a/4/ (27/3) — l)/3. In the laboratory frame, these 



two solutions describe sinusoidal oscillations with fre- 
quency ujp, dimensionless oscillation amplitude \Q\\ and 
\Qh\, with |Qi| < \Qh\, and opposite phases, <f)\ = and 
4>h = ft- They describe the low- and high-amplitude so- 
lutions, respectively. 

For weak damping, k = Y / \8ui\ -C 1, the low- and high- 
amplitude states become attractors. For weak quantum 
fluctuations, A« 1, there is a separation of time scales: 
On the timescalc ~ re -1 , the oscillator relaxes to the 
vicinity of either one of the two states, while on a much 
longer timescale oc exp[A _1 ] rare large thermal or quan- 
tum fluctuations induce switching between them. In the 
next section, we focus on the quasi-stationary regime, 
k _1 <C r <^ K _1 cxp[A _1 ], where the oscillator stays in 
the vicinity of one of the two metastable solutions. 



B. Small fluctuations around the mean field 
solutions 

In the weak damping regime, k<1, small quantum 
and thermal fluctuations around a mean field solution can 
be taken into account in two steps. In the first one intro- 
duces an auxiliary oscillator, which describes the quan- 
tum states localized in the phase-space region close to 
one of the two solutions. In the second one introduces an 
effective master equation for the auxiliary oscillator. It 
turns out that the fluctuations and the dissipation for the 
auxiliary oscillator are similar as for an oscillator close to 
equilibrium but with an effective temperature T e , to be 
discussed below Q. 

The auxiliary oscillator is obtained by expanding the 
quasi-energy Hamiltonian © around the classical solu- 
tions <3; = Qi/h and P = 0, where the index i denotes 
the states of low- (i=l) and high- (i=h) amplitude oscil- 
lations: 



9 ~ 9i + -j9ppP 2 + ^9qq{Q - 



(11) 



with g { = g(Qi, 0), g PP = Q? - I and g QQ = 3Q? - I. 

It is convenient to rewrite the Hamiltonian in terms of 
ladder operators 

9 « 9(Qi, 0) + sgn( 3QQ )A v x {b% + 1/2) (12) 

where v x = JgQqgpp is the dimensionless frequency of 
the slow oscillations and the ladder operators b and 6* 
are obtained by the squeezing transformation J19| 



CJi 



b cosh rt — fo T sinh i 



(13) 



with ai = Qi/v2A and squeezing factor r* = 

M5qq/3pp]/4- 

For k< Co we can incorporate small thermal fluctu- 
ations by substituting the squeezing transformation Eq. 
(fTB")) into the Lindblad master equation © and neglecting 
the terms where the operators b and b' are not matched 
in pairs. We find 

d T p = -j\g,p] + k(1 + n e )V[b]p + K n e T>[tf]p, (14) 



with n e = n+ (2n + 1) sinh 2 r* [5[. This master equation 
has the same structure as the master equation for a static 
weakly-nonlinear oscillator with eigenfrequency v\ and 
damping k at the temperature T c = \v\/\u[{n e + l)/n e ]- 
It naturally yields a stationary distribution for the aux- 
iliary oscillator Fock states In) of the Boltzmann form 

SGI 



p st = 

rn — 



{n\p st \n} pa (n e + l)" 1 exp[-nA^/T e ] 



(15) 



Note that for a driven oscillator, the effective tempera- 
ture T c is finite even for n = 0, T c = Afi/(21ncothr*) [5j, 
Il4 , |2CJ , aphenomenon which has been named quantum 
heating Q- In the classical limit of high temperatures 
ksT ^> hajQ, we have [5| 



T e « A^in e K,k B T 



■il 



Auj^Suj 



(l + 2sinh 2 r*), 



(16) 



so that T e becomes independent of H, as it should be 
expected for a classical quantity. 



IV. THE EMISSION SPECTRUM 

The master equation for the auxiliary oscillator Eq. 
([T4"]) allows us to compute the emission spectrum S(u) 
close to the Lorentzian peaks, in the limit where they are 
well separated, k< vj. In general, it is not sufficient to 
use the harmonic oscillator approximation for the auxil- 
iary oscillator, rather one has to include in Eq. ([T4"|) the 
relevant nonlinear corrections to g. To proceed, we focus 
on the parameter regime where the second order peaks 
are clearly visible. In this case, we can neglect in the 
quasi-energy Hamiltonian g (J5J terms which are quartic 
in the ladder operators b and w . Thus we find 

g « $ + sgn{g QQ )\ Vi tfb + A 3 / 2 (^6 3 + V 2 tf 2 b + h.c.) , 

(17) 
with Vi = Qi(e- 3r ' -e r *)/2 3 / 2 and V 2 = Qi(3e- 3r . + 
e<)/2 3 / 2 . 

In order to compute the emission spectrum S(w), we 
insert the squeezing transformation Eq. (|13l) into the def- 
inition Q of the emission spectrum and obtain the ex- 
pression 

S(u) = TT^-r (cosh 2 rf 2Re(6 f , b)„ + sinh 2 r*2Rc(6, tf) v 
\5u\ 



— 2sinhr* coshr*2Re(6, b) u ) . 
Here, we have defined 

v = (u> — u)f)/\5u}\ 
and the spectral functions 

(0,0')„= / dre~ WT (0(r)0') st 



(18) 



(19) 



(20) 



The evaluation of the auxiliary oscillator spectra (b, b') v , 
(b\b) v and (b,b) v proceeds along the same lines as the 



calculation of the spectrum of a weakly-nonlinear oscilla- 
tor [l3:[2l|, and is detailed in Appendix 1X1 By substitut- 
ing Eqs. (|A10j> . (|A11|I and (|A12|) into Eq. (JTSJ we obtain 
the results for S(w) close to its peak. 

For frequencies v close to v\ and — t'; we find the first 
order peaks 

2k cosh r*n e 



and 



S(w) = 



S(w) = 



f 
|<5u;| (1/ - sgn(.gQQ)i/i) 2 + re 2 ' 

f 2Ksinh 2 r*(n e + l) 



(21) 



(22) 



\6w\ iy + sgn( ff Q Q )i/i) 2 + k 2 ' 

respectively. The expressions for the first order peaks 
have been obtained previously using a quantum Fokker- 
Planck equation approach [la ] . 

For frequencies v close to zero we find 



5(w) 



|<5w| v 2 { ' + J ^ 2 +4k 2 



(23) 



The central peak is a consequence of amplitude fluctu- 
ations Sa n = (n\a\n) — (a) st ^ of the quasi-energy 
states. It is in fact possible to rewrite eq. ([2"3"[) as 
2(d5a n /dn) 2 {tfb,tfb) u [HJ. The peak is therefore de- 
termined by a combination of the amplitude curvature 
and quasi-energy fluctuations. Close to ±2i/o we find the 
second-order peaks 



5(w) 



S(u) = 



1 XQf 



16k cosh 4 r* n 2 



\Su\ v?e Ar ? (v - sgn(<7Q Q )2^) 2 + 4k 2 
1 XQ 2 16Ksinh 4 r*(fi e + l) 2 



(24) 



(25) 



\Su\ v?e Ar ! {v + sgn(<?Q Q )2^) 2 + 4k 2 

Equations ([2"5[) . (J2U) and ([231) are valid only close to the 
corresponding peaks, which do not have a substantial 
overlap for k 2 <C A. Moreover they do not incorporate 
a shift in the peak center of order A|(5u;|, and do not 
take into account the non equidistance of the quasicncrgy 
levels, which give rise to a hyperfine splitting similar to 
Refs. \M and U3, when it exceeds the levels broadening, 
for k <C A. Therefore strictly spcking they are valid for 
k 2 < A < k. 

In the classical limit fc^T ^> hajo, the intensity of the 
emitted power irradiated ThuS(u}) becomes independent 
of h as can be verified by substituting Eqs. (|5I16[) into 
the lineshape Eqs. ([2T1 - |2"5[) of S(oj) close to its peaks. 

The amplitude of the first order peaks is linear in the 
quantum temperature h e , while the second order peak 
are proportional to the squared noise intensity n 2 . Fig. 
[3] shows the quantum temperature n e of the two sta- 
ble states for n = 0. As we can see, in the state of 
low-amplitude oscillations, the quantum temperature is 
mostly rather small. It becomes larger for high values 
of the scaled driving strength /3; however, in this regime, 
the stationary state of the driven resonator is the state of 
high-amplitude oscillations. On the other hand, for high- 
amplitude oscillations, the quantum temperature takes 
large values. 
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Figure 3. The quantum temperature n e for the two stable 
states as function of the scaled driving strength /3 for n = 
(see Ref. [3). 



V. RESULTS 

In this section we present the results for the emis- 
sion power spectrum. We mostly show results from a 
numerical solution of Eq. ([3]) with boundary condition 
p(0) = ap s , which also covers parameters beyond the 
regime k 2 <C A <C k, where the analytic solution is valid. 

First we focus on the case of strong driving, where 
the resonator is in a stationary state of high-amplitude 
oscillations. In Fig. [4] we show the emission spectrum for 
finite temperature, n — 1.0 and /3 = 0.13. In the regime 
of high-amplitude oscillations, the squeezing factor r£ is 
rather large, so that according to Eqs. (|21[) and (|22l) . 
both first-order peaks are of comparable size. 

We observe the appearance of two second-order peaks. 
The first one is located at the frequency u> w — 2 \6u\ 
and is described by Eq. ([2~4")l . The oscillation amplitude 
Qh is large, \Q^\ m 1, so that the peak height ~ X/k is 
larger than the background ~ k arising from the first- 
order peaks. We also observe a second-order peak at the 
frequency ui = 0, which is described by Eq. ([2U[) . We note 
that for A = 0.005, where the condition k 2 < A < k is 
fulfilled, the numerical results are very well described by 
the analytical expressions ([2"Tl) - ([2"5"[) . 

The analytical expressions we derived for the emis- 
sion spectrum predict that the height of the second-order 
peaks is proportional to the ratio X/k. Therefore the 
second-order peaks become more visible against the back- 
ground for larger values of the parameter A. This is also 
shown in Fig. 31 where we plot the emission spectrum for 
different values of the parameter A. 

We see that, while the first-order-peaks are decreasing, 
both second-order peaks are clearly increasing for larger 
values of A, consistent with the analytical predictions. For 
the large values of A used in the plot, we are close to the 
limit of validity of the analytical expressions for the spec- 
trum, since the condition A <C k is not strictly fulfilled. 
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Figure 4. (a) Emission spectrum of the driven nonlinear res- 
onator in the regime of low-amplitude oscillations for three 
different values of the parameter A. The other parameters are 
/3 = 0.13, k — 0.02 and n = 1. (b) and (c) shows a zoom 
around the second order peaks at the frequencies ui « — 2 \8ui\ 
and u) — 0. 



In this regime our analytical result still provide the cor- 
rect qualitative picture for the emission spectrum and a 
remarkably good estimate of the line intensities, but it 
does not describe a finite shift of the peaks positions. 

For the sake of comparison, we also discuss the case of 
weak driving , where the resonator is in a state of low- 
amplitude oscillations. The emission spectrum is shown 
in Fig. [5] for zero and finite temperature, and for fixed 
driving, /3 = 0.01. Here only the first-order peaks, lo- 
cated at the frequencies u ps ± \6u>\ are visible. 

Since the oscillation amplitude Q\ is small, Q\ -C lj 
the squeezing factor r* is close to zero, r* « 0, and the 
effective temperature is close to the physical tempera- 
ture, n e « n. This means that for finite temperature, 
n > 0, the first-order peak ~ n e cosh r* at the frequency 
uj w | (5a; | is much larger than the second first-order peak 
~ (n e + 1) sinh r* at the frequency u> ss |(5o;|. 

On the other hand, in the limit of zero temperature, 
when n e tends to zero, both first-order peaks are of 
the same (small) size. Finally, the second-order peaks 
~ XQq are so small compared to the background of the 
first-order peaks that they are not visible. 
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Figure 5. Emission spectrum of the driven nonlinear resonator 
in the regime of low-amplitude oscillations. The parameters 
are /3 = 0.01, A = 0.005 and k = 0.02. 



VI. SUMMARY 

We have studied the emission spectrum of a driven 
nonlinear oscillator. We focused on a regime character- 
ized by a strong nonlinearity and low temperature, which 
are typical for circuit QED experiments and lead to the 
appearance of pronounced nonlinear effects. 

In the case that the oscillator is locked to the regime 
of high-amplitude oscillations, we found two new fea- 
tures in the spectrum, in addition to the Raman-side 
peaks located at frequencies u ~ ± |<5u;|: i) We observed 
the appearance of a second-order peak at the frequency 
u> ~ u>f — 2\5u\. In a quantum-mechanical picture, this 
emission peak results from transitions between next-to 
nearest quasi-energy levels, ii) We observed a second- 
order peak at the driving frequency up (in addition to 
the sharp Rayleigh peak) which results from fluctuations 
in the asymmetric quasi-energy potential. 

Based on a master equation approach, we derived an- 
alytical expressions for the emission spectrum around 
the first- and second- order emission peaks. The results 
were obtained by expanding the quasi-energy Hamilto- 
nian around the classical stable solutions. In order to de- 
scribe the second-order peaks, we went beyond the usual 
linearization of the equations of motion and accounted for 
the nonlinear nature of the quasi-energy Hamiltonian. 
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Appendix A: Calculation of the spectra for the 
auxiliary nonlinear oscillator 

In this Appendix we compute the spectra (b, tf) v , 
(b\b) u and (b,b) v , for the auxiliary nonlinear oscillator, 
defined in Eq. (fT7)) . whose dissipative dynamics is gov- 
erned by the Lindblad master equation (|I4|) . We start 
with 

/>oo 

(&,&*>„= / dre- iVT {b{T)tf) st (Al) 

Jo 

By integrating by parts we get 

(b,tf) v = --(b,tf) v -~0). 



(A2) 



The equation of motion for (b), (b) = i/X([g,b\) — k(&) 
[19J . directly applies to (b(t)b) s t, since (b(t)b) st = (b(t)} 
with initial condition (b(0)) = ti(b^bp st ). By inserting it, 
together with the approximate expression of the quasi- 
energy g into Eq. (fl7)) , into Eq. (|A2[) . we find 



(M f >, 



{B,tf) v + i{btf) 
v + sgn(g QQ )v - %k 



(A3) 



with B = A^pVifet 2 + 2V 2 b% + V 2 b 2 }. Taking into ac- 
count that 

{B,b) v = dTe-™ T (B(T)tf) st = dre™ T (Bb\T)) st 

JO J -oo 

(A4) 
we can integrate by parts and apply the equation of mo- 
tion for (b^), obtaining 

(,,0,-^m. (A 5 ) 

z/ + sgn(g Q QJ^o + «'>: 
By substituting into Eq. (jA3|) . we arrive to the identity 
t = (B.Bt^ + ifBfct) j(w,t) 



(z/ + sgn(.gQQ)^ ) 2 + k 2 ^ + sgn(g QQ )vo - in 

(A6) 
Next we replace the nonlinear spectra in the right hand 
side of Eq. (|A6|) with the spectra of a weakly-damped 



harmonic oscillator 
(b 2 MX 
(rf 2 ,b 2 ) ^ 



AK{fi e + l) 2 



v + sgn(gQQ)2v - i2n 

4kti 2 



V 



sgn(ffQQ)2i/ - i2n ' 



f — i2n 



(A7) 
(A8) 
(A9) 



and find 



(bM)u 



n(n e + 1) 



v + sgn(gQQ)^ - i« vl 



XVi (b\brx 



-^\bX + ^{b\b%)l. (AIO) 

"0. ^0 

The result consists of four terms, each of them yielding 
a Lorentzian peak. It is important to keep in mind that 
the substitution of a nonlinear spectrum with a linear one 
is valid only close to the peak of the spectrum. Hence, 
our approximate result for (b,b^) v gives a good approx- 
imation close to its second order peaks at ±2i/o,0, only 
when the peaks are well pronounced. The peak height, 
given by Eq. (|A10[) . is of the order A/n. On the other 
hand the background value can be computed by lineariz- 
ing the fluctuations, and is found to be of order k \W^ . 
This leaves us with the requirement k 2 <C A. 

The calculation of the spectra {b\b) v and {b,b) v pro- 
ceeds along the same lines and yields 



{b\b) v 



^t(b\bV)l 



v - sm(g QQ )v - in v^ 
+M (6 t 2;62) o + A4^ (6 t 6i6 t 6) o. (An) 



and 



(6,fe)^-^((& 2 ,^)0 + (& t 2j6 2 ) o ) 
v o 

+^L (6 t 6j6 t 6> o. 



(A12) 
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